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Abstract: Parton shower Monte Carlo event generators in which the shower evolves from 
hard splittings to soft splittings generally use the leading color (LC) approximation, which 
is the leading term in an expansion in powers of 1/A^, where Ac = 3 is the number of col¬ 
ors. In the parton shower event generator Deductor, we have introduced a more general 
approximation, the LC-|- approximation, that includes some of the color suppressed contri¬ 
butions. In this paper, we explore the differences in results between the LC approximation 
and the LC-|- approximation. Numerical comparisons suggest that, for simple observables, 
the LC approximation is quite accurate. We also argue that for gap-between-jets cross 
sections neither the LC approximation nor the LC-|- approximation is adequate. 
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1 Introduction 

Parton shower Monte Carlo event generators like Pythia [1], Herwig [2], and Sherpa 
[3] represent the momenta and flavors carried by partons in a natural way [4-7]. One can 
imagine that at any “shower time” t in the development of the shower, there is a probability 
p{t,{pa, fa,ph, fh,Pi, fi, ■ ■ ■ ,Pm, fm}) to have 771 flnal-state partons plus two initial-state 
partons with the specifled momenta pi and flavors fi. What happens between time t and 
t -b At is a change A\p{t)) chosen at random according to physics based parton splitting 
probabilities. The evolution of |/?(t)) is straightforwardly described in the language of 
classical statistical mechanics. 

Spin and color do not fit easily into the event generator format because quantum 
interference between different spin and color states is important. For this reason, the 
natural language is that of quantum statistical mechanics. Then |/o(t)) represents not just 
a probability distribution, but a density matrix in color and spin space. 

There is a fairly straightforward formalism [8] available to specify an evolution equation 
for |/?(f)) that takes full account of color and spin. Unfortunately, we do not know how 
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to implement the full evolution equation in a practical computer program that would 
constitute a parton shower event generator. If we use the leading color (LC) approximation 
and average over spins, we do get something practical [9]. 

We have implemented the resulting shower evolution equation as a parton shower event 
generator, Deductor [10, 11]. In ref. [10], we presented results from Deductor using the 
leading color approximation. However, there is an improved version of the leading color 
approximation that is still practical: the LC+ approximation [12]. The LC+ approximation 
is built into Deductor. 

In this paper, we use the LC+ approximation to explore numerically the effect of color 
in a parton shower. We do not have a shower treatment with full color. It should be 
possible [12] to expand cross sections perturbatively in the difference between full color 
splitting functions and their LC+ approximations, but this possibility is not implemented 
in Deductor. Thus we are limited to looking at what we can discover with the LC+ 
approximation. From this investigation, two tentative conclusions emerge. 

First, comparing results for simple observables calculated with the LC and LC+ ap¬ 
proximations, we hnd differences of order l/Al^ ~ 10%, but not larger differences. This is 
perhaps not a surprise, but real calculations were needed to reach this conclusion. 

Second, in events with two high px jets with a large rapidity separation between them, 
we argue that there can be a substantial influence of color on the probability to have a gap 
between the two jets that contains no jets with transverse momentum above a cut p™*. 
Here, neither the LC approximation nor the LC-|- approximation is adequate. 

This paper is organized as follows. We provide some background information in section 
2. Then in section 3 we investigate how much color suppression is typical in the color state 
as the shower develops. (The brief answer is that it is rare to have no color suppression.) 
In sections 4 and 5, we look at the difference between LC results and LC-|- results for two 
simple observables: the one jet inclusive cross section and the distribution of the number 
of partons in a jet (which is not really physically observable, but is calculable with a fixed 
cutoff on the /ct in parton splittings). In section 6, we examine the influence of color on the 
cross section to have a gap between high px jets. We present some conclusions in section 
7. There is an appendix A about how to switch from the LC-|- approximation to the LC 
approximation part way into the shower. Another approach to color in a parton shower 
has been presented in the literature [13] . We compare this approach to ours in appendix B 

2 Background 

Before exploring numerical results, we provide a little background about Deductor and 
about color and the LC-|- approximation. 

2.1 Deductor 

A brief presentation of the parton shower event generator Deductor can be found in 
ref. [10]. Here, we simply mention some major features. First, Deductor generates a 
parton shower, but does not at present include a hadronization stage. Nor does it include 
an underlying event. Thus it is well suited as a tool for investigating the approximations 
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in parton shower algorithms, but has some limitations for other uses. Second, Deductor 
averages over parton spins at each stage of splittings, in the same fashion as other parton 
shower event generators. There is a method for including spin that should be practical 
[14], but we have not yet implemented it. Third, we argue in ref. [15] that it is best to 
order splittings in a parton shower in order of decreasing values of the virtuality in the 
splitting divided by the energy of the mother parton. This rather non-standard choice 
is implemented in Deductor. Fourth, initial-state charm and bottom quarks have their 
proper masses in Deductor [16]. 

The parton splittings in Deductor are ordered from hardest to softest. Quantum 
interference is accounted for at least approximately because the shower is based on color 
dipoles [17, 18]. Another approach is to order splittings in angles, from largest to smallest 
[6, 7]. In this approach, one makes a rather severe approximation with respect to the 
dependence of the splitting probability on the azimuthal angle of the emission, but one 
gains accuracy with respect to color structure. This is the approach followed in Herwig 
[ 2 ]. 

2.2 The color density matrix 

The natural language for describing an evolving probability distribution in a parton shower 
is statistical mechanics. In order to include quantum color, we need quantum statistical 
mechanics. Thus we need a density operator that depends on the shower time t and is an 
operator on the space of color states of a possibly large number of partons. The density 
operator has the form 

P{{pJ}m,t)= ^ P{{P,f,c',c}m,t)\{c}m){{c'}m\ ■ (2.1) 

Here |{c}m) and |{c'}m) are standard [19] basis vectors for the quantum color space for m 
final state partons plus two initial-state partons [8, 12]. The color configuration {c}m of the 
ket state is, in general, different from the color configuration {c'}m of the bra state. Thus 
the function p{{p, f,c' depends on two sets of color indices. The density operator 
p{{p, can be regarded as a vector in the space of functions of {p, f}m with values 

in the space of operators on the quantum color space. 

The color basis states are normalized to {{c}m\{c}m) = 1 or to {{c}m\{c}m) ~ 1 with 
very small corrections. They are not, however, generally orthogonal. However, when {c}m 
and {c'}m are different, one finds that {{c'}m\{c}m) = with P > 1. That is, the 

basis vectors are orthogonal in the Nc ^ oo limit. 

At the end of a parton shower, one will want to measure something about the final 
state, using a color singlet measurement operator.^ This means that we take the trace 
of p{{p, f}m,t) in the color space, so that the contribution from color states {c',c}m is 
proportional to {{c'}m\{c}m)- That is, the contribution from color states with {d}m 7 ^ 
{c}m is suppressed by one or more factors of 1/iVc. 

^This leaves out hadronization. If we apply a hadronization model based on the formation of color 
strings, we do, in effect, apply a non-color-singlet measurement operator. We formulate how this can work 
in ref. [12]. 
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One notes, however, that states with {c'}m / {c}m are a natural part of quantum 
chromodynamics (QCD). They appear even at the Born level of a hard scattering, where 
we normally keep track of them by using color ordered amplitudes. Even if we start with 
{c'}m = {c}m states, splittings in a parton shower generate {c'}m / {c}m states. 


2.3 The color suppression index 

In order to understand the systematics how many factors of 1/A^c accompany states with 
{c'}m 7^ {c}ni, it is useful to use the color suppression index [12]. 

To do that, we first note a small technical complication. In a g 
have a color matrix in the color amplitude. We can use the Fierz identity. 


q + q splitting we 
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( 2 . 2 ) 
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to write the result expanded in color basis states. At each g ^ q + q splitting, shower 
evolution picks either the first, leading color, term or else the second, color suppressed, 
term. If the second term is chosen, further evolution uses the second color state, 
and incorporates the factor — 1/(2W) into the weight factor for the event. We let pE 
represent the number of times during the shower evolution that we pick up a l/W factor 
by using the second term in the Fierz identity. 

With that complication out of the way, we can define the color suppression power P as 
the number of powers of 1/Nc that appear at some stage of the shower coming from both 
the overlap ({c'}m|{c}m) of the color states, and the explicit 1/W factors from g ^ q + q 
splittings: 


/ 1 \pe 
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Now we define the color suppression index I. We simply use the group U (W 
of the true color group SU (W) to calculate the color overlap 
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(2.3) 
in place 

(2.4) 


The color suppression index has two properties that make it quite useful. First, we 
always have 

P{m) > I{m) . (2.5) 

That is, if we use I{m) to estimate the amount of color suppression, we can never overes¬ 
timate. Second, at each stage of the shower, the color suppression index either stays the 
same or else it increases: 

I{m + 1) > I{m) . (2.6) 

Thus we can think of I as measuring color disorder, like entropy: it can never decrease as 
the shower evolves. 


2.4 The LC and LC+ approximations 

We have noted above that the hard scattering that leads to a parton shower has some 
contributions with {d^m = {c}m and thus with color suppression index I = 0.^ It also 

the hard scattering is 2 —>■ 2 parton-parton scattering, then m = 2 at the start of the parton shower. 
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generally has some contributions with {c'}m 7^ {c}m and thus with color suppression index 
7 > 0. This implies that the associated probability has at least I powers of 1/iVc. As the 
shower develops, each parton splitting gives some terms that leave I unchanged, but also 
some terms that increase I. Thus, inevitably, we generate some color conhgurations that 
come with high powers of 1/iVc and thus low probabilities. 

The contributions with 7 > 0 are also associated with calculational difficulties in a 
parton shower. For this reason, it is common to start at the hard scattering with {c'}m = 
{c}m, 7 = 0 states only and then to omit generating any 7 > 0 states. In practice, it is 
very easy to avoid increasing 7 in parton splittings. In simplest terms, we count each gluon 
as carrying color 3x3 instead of color 8. This gives us the leading color approximation. 

We can also use the LC+ approximation [12]. Here we can start with {c'}m /{c} mi 
7 > 0 states from the hard process. Then we can keep some of the terms with increasing 
7 that are generated by parton splittings. We do not keep all of the terms, so this is still 
an approximation. The approximation is exact for collinear splittings and it is exact for 
splittings that are simultaneously collinear and soft. It does not, however, keep everything 
that arises from wide angle soft splittings. We refer the reader to ref. [12] for a full 
description. 

With the LC+ approximation, we can eventually generate states with quite high values 
of the color suppression index 7. Because these states come along with numerical weights, 
this has the effect of slowing down the numerical convergence of the program. That is, it 
takes more events to give the same statistical error. For this reason, we set a maximum 
value (A7)max for how much can be added to the color suppression index by the parton 
shower. When the amount of color suppression added by the shower, A7, reaches (A7)inax, 
we switch off the LC+ approximation and start using a version of the LC approximation 
that is adapted to {c'}m / {c}m- We describe how to do this in appendix A. For numerical 
results in this paper, we set (A7)max = 4. Shower splitting either leaves 7 unchanged or 
increases it by two. Thus we are omitting terms with A7 = 6, corresponding to omitted 
factors as large as 1/N^ ~ 10“^, which should allow for sufficient accuracy for any practical 
purpose. Of course, (A7)niax is adjustable, so that we have an approximation with errors 
that can be controlled. 

2.5 Systematic improvement in color 

Let us look briefly at how the LC+ approximation works and how it can be improved [12]. 
With full color, shower evolution has the form 

\p{t)) =U{t,to)\p{tQ)) , (2.7) 

where \p{t)) represents the state of the system at shower time t in a statistical ensemble 
of trials. The evolution operator 7/(t,to) obeys the evolution equation, 

^ U{t, to) = [ni{t) - V(t)] U{t, to) . (2.8) 

Here 'Hi{t) generates a parton splitting, while V{t) leaves unchanged the number of partons 
and their momenta. For the LC+ approximation, we approximate U{t,to) by {t,to) 
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where 


(2.9) 


This differential equation can be solved iteratively in the form 

U^^+{t,to) = M^^+{t,to) + [ dTU^^+{t,T)'H^^+{T)Af^^+{T,to) . ( 2 . 10 ) 

Jto 

Here is the no-splitting operator, 


AA^'^+(t 2 ,ti) = exp 



dr V^‘^+(t) 


( 2 . 11 ) 


That is, Af^^~^{t 2 ,ti) is the Sudakov factor that represents the probability that no parton 
splits between shower time ti and shower time t 2 - The essential point of the LC-|- approx¬ 
imation is that the operator V^'^'’'(r) is diagonal in the standard color basis that we use, 
so that it is practical to calculate its exponential. In contrast, the operator V(r) is not 
diagonal in color, so that it is not practical to evaluate its exponential. 

Now, what if we want shower evolution with full color? Then we need 


+ Ani{t) - Avit)]uit,to) , ( 2 . 12 ) 

where 

Ani{t) = niit)-n^^+{t) , 

AVit) = Vit)-V^^+{t) . ^ ^ ^ 


This evolution equation gives 

l({t,to) + f dr ??(t,r) [A?^/(t) - AV(r)]??^*^+(r,to) , (2.14) 

Jto 

which can be solved iteratively: 

U{tuto)=U^^+{ti,to) 

+ f dr [Aniin) - AV(ri)] ??^'^+(ri, to) 

Jto 

+ f dT 2 f dn l{^^^{ti,T 2 ) [Ani{T 2 ) - AV{t2)]1(^^^{t2,Ti) 

J to J to 

X [A?t/(Ti) - AV(ri)] U^^^{Ti,to) 

+ ■■■ . 


One can have any (small) maximum number of insertions of [A'H 7 (r) — AV(r)]. In between 
the insertions, we have LC-|- evolution. The operator [A??/ (r) - AV(r)] can change the 
color state of the system, but we can still proceed because {TiJ^i^Ti) can act on any 
color state. 

We will argue in section 6 that this sort of systematic improvement is needed for 
some problems. We have, however, not implemented this feature in the current version of 
Deductor. 
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3 How much color suppression is typical? 

Is it rare to generate a state with non-zero color suppression index /? To find out, we 
use Deductor to generate the cross section to make jets in proton-proton collisions at 
^/s = 14 TeV. We select events that have a jet with pT > 500 GeV in the rapidity range 
-2<y<2. The parton shower generates splittings that are softer and softer, as measured 
by the virtuality parameter defined in ref. [15]. We stop the shower by not allowing any 
splittings with a splitting transverse momentum smaller than /dp™ = 1 GeV. This allows 
for quite a number of splittings: most of the cross section comes from events with twenty 
to fifty splittings. We show the distribution of the number of splittings in figure 1. 

The hard process here is parton-parton scattering.^ We keep all of the contributions 
to the color-ordered amplitudes for parton-parton scattering. This gives an initial color 
density matrix of the form (2.1) with m = 2 for the two final state partons. We calculate 
the color suppression index I2 just after the hard scattering according to eq. (2.4) with 
Pe = 0. We find that, including the cases with {c'}2 7^ {c}2, the color suppression index 
can be I 2 = 0, 1, or 2. Now the shower starts. With each splitting, I has a chance to 
increase, so that Im > h- We set a limit Im — h ^ (A/)max with (Al)max = 4:. If Im — h 
reaches /max, we switch off the LC-|- approximation and start using the extended version 
of the LC approximation that allows for {c'}m m* 

We can answer the question of how rare it is to generate a non-zero color suppression 
index by plotting the fraction of the jet cross section coming from events that have a given 
value of In+2 after N splittings. We group In+2 = 1 with In+2 = 2, In+2 = 3 with 
In+2 = 4, and In+2 = 5 with In+2 = 6 because odd values of I (which are generated at 
the hard scattering) are rather rare. We plot the fraction of events with given values of 
In+2 for N up to 50. For events with S splittings with S < 50, we define In+2 = Is+2 
for N > S. We show a plot of the fraction of events with given values of /Ar +2 after N 
splittings in figure 2. 

We can draw two lessons from figure 2. First, a little less than 10% of the cross section 
is associated with events that start with / > 0 right at the hard scattering. This is about 
what one would expect, since 1/N^ ~ 0.1. Second, it is rather rare for the partons to stay 
in an / = 0 state after many splittings. Only about 20% of the cross section comes from 
parton configurations that have / = 0 at the end of the shower. In fact, more than half 
of the cross section has / > 2 at the end of the shower. We would see more cross section 
associated with / > 4 at the end of the shower if we had not imposed /^ — /2 < 4 on the 
shower. 

4 Jet cross section 

We have seen that most of the cross section to make jets in proton-proton collisions at 
^/s = 14 TeV is associated with partonic states at the end of the parton shower that have 
color suppression index / > 0. Does this matter for the jet production cross section? 

®We have checked our calculation of the color density matrix for parton-parton scattering against the 
results of Kuijf [20]. 
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Figure 1. Probability of generating N splittings versus N in events that have a jet with px > 
500 GeV in the rapidity range —2 < y < 2. The splitting cutoff is fc™™ = 1 GeV. 

One would expect that it does not matter, since the jet production cross section is so 
inclusive. That is, by design, the cross section to make a jet with a certain px in a rapidity 
range —2 < y < 2 is very much infrared safe. The cross section is mostly determined by 
the cross section for parton-parton scattering (calculated at the Born level in this case). 

Barton shower evolution can add some px to the jet from partons radiated from the 
initial-state partons. Shower evolution can also subtract some px from the jet when the 
scattered parton that initiates the jet radiates a gluon at such a wide angle that the 
gluon is not included in the jet. This initial-state and final-state radiation, along with 
the interference between initial-state and final-state radiation, is influenced by the color 
configuration of the radiating partons. Thus the development of the color configuration 
can influence the ultimate jet cross section. Since the color state develops differently in 
the LC-b approximation than it does in the LC approximation, the jet cross section could 
be affected. However, since the jet cross section is designed to be insensitive to all of these 
infrared or long-time effects, one does not expect to see much difference between the jet 
cross section calculated in the LC-|- approximation and the jet cross section calculated in 
the LC approximation. 

We can test the hypothesis that the LC approximation gives a jet cross section close 
to that of the LC-b approximation. We use Deductor to calculate the one jet inclusive 
cross section da/dpT to produce a jet with transverse momentum px in the rapidity range 
—2 < y < 2. With the help of Fast Jet [21], we use the /cx jet algorithm [22, 23] with 
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Figure 2. Fraction of events with given values, /, of the color suppression index as a function of 
the number of splittings. We set a limit of 4 on the change of I after the hard scattering. 


R = 0.4. The results are shown in figure 3. The solid line shows the cross section calculated 
with the LC+ approximation with (A/)inax = 4. The dashed line shows the result obtained 
with the LC approximation.^ The dashed line and the solid line are indistinguishable on 
the graph. We would show a graph of the ratio of the two cross sections, but we find that 
the difference is less than 2%, which is the statistical accuracy of the calculation. 

We conclude that the use of the LC approximation for the shower instead of the LC+ 
approximation makes a negligible difference for the one jet inclusive cross section, even 
though most of the cross section comes from states that have color suppression index / > 0 
by the end of the shower. This is not a surprising conclusion, but it is not a conclusion 
that one could be sure of without doing the calculation. 

5 Number of partons in a jet 

In section 4 we looked at an observable that is very insensitive to soft parton splittings: the 
one jet inclusive cross section. Now, we look inside these jets at a quantity that is sensitive 
to soft parton splittings: the distribution of the number of partons in a jet. Evidently, the 

^In the hard scattering, density matrix contributions |{c} 2 )({c'} 2 | with {c }2 7 ^ {d }2 are eliminated 
and the corresponding cross section is redistributed to density matrix contribntions |{c} 2 )({c'} 2 | with 
{c }2 = {c'} 2 - Then (A7)max is set to zero in the shower so that the color snppression index stays eqnal to 
zero. 
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Figure 3. One jet inclusive cross section da/dpT for \y\ < 2 using the fcx algorithm with R = 0.4. 
The red solid line shows the cross section calculated with the LC+ approximation. The blue dashed 
line shows the cross section calculated with the LC approximation for the parton shower. The two 
curves are indistinguishable. 

number of partons in a jet is not a physical observable, but it is calculable with a fixed cut 
on the smallest transverse momentum allowed in a parton splitting. It is a stand-in for the 
number of hadrons in a jet, which is physically observable. The distribution of the number 
of partons in a jet is of interest here because it is an infrared sensitive quantity. 

We analyze a sample of jets with pp > 200 GeV and \y\ < 2. We examine the distri¬ 
bution Pn{n) of the number n of partons in a jet in this sample for events simulated by 
Deductor using the LC-|- approximation and using the LC approximation. The distribu¬ 
tion is normalized to Yin ~ each case, we stop the shower by not allowing any 

splittings with a splitting transverse momentum smaller than /c™™ = 1 GeV. 

The result is displayed in figure 4. We see that for n < 10, Pn{n) is quite insensitive 
to whether we use the LC or LC-|- approximation. The behavior of Pnin) for large n is 
of special interest because when there are many partons, most of them must have small 
momenta, so that for large n we probe the softest emissions. The probability that there 
are more than 10 partons in a jet is small. We notice that for n > 10, Pn{n) calculated 
using the LC-|- approximation is about 10% larger than pn{n) calculated using the LC 
approximation. The fractional difference is of order 1/V^, in accordance with the simplest 
expectation. 
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Figure 4. The distribution pn{n) for the number n of partons in a jet with > 200 GeV and 
\y\ < 2. We compare the distribution calculated with the LC+ approximation (solid red histogram) 
with the with the same distribution calculated with the LC approximation (dashed blue histogram). 
The results are very close for n < 10. The jets are constructed using the A:t algorithm with R = 0.4. 

6 Gap fraction 

We now turn to an observable for which we believe that the color treatment matters: the 
probability that there are no jets in the rapidity range between two jets that are widely 
separated in rapidity. 

In proton-proton collisions at ^/s = 14 TeV, we look for jets using the algorithm 
with R = 0.4. We select events for which there is at least one jet with transverse momentum 
Pt > 200 GeV and rapidity y > 2 and at least one jet with transverse momentum px > 
200 GeV and rapidity y < —2. For a given transverse momentum p™*, we can ask for the 
fraction /(p™*) of selected events in which no jet that has transverse momentum greater 
than p™* appears in the rapidity region —2<y<2. 

There has been considerable work on the theory of this sort of observable [24-34] . The 
treatment of the problem in an angle ordered parton shower is analyzed in ref. [35]. 

In figure 5, we plot /(plp*), versus p™*. For large values of p™*, the probability to 
produce an extra jet is suppressed by a factor ccg and by the need to have parton distribution 
functions with larger momentum fraction arguments, so that we expect /(p™*) to be close 
to 1. However, it is typically rather easy to produce low transverse momentum jets, so 
for small p™'', we expect /(p™'') to be small. This is what we see in figure 5. We have 
calculated /(p™*) both in the LC approximation and in the LC-t- approximation. We see 
that there is about a 10% difference between the two results at small values of p™*. This 
level of difference is similar to what we saw for the distribution of the number of partons 
in a jet in section 5. 
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Figure 5. For events with at least one jet with transverse momentum > 200 GeV and rapidity 
y > 2 and at least one jet with transverse momentum px > 200 GeV and rapidity y < —2, we 
plot the fraction of events that have no jets with transverse momentum greater than p™* 

in the rapidity range —2 < y < 2. The blue dashed curve shows the result obtained with the LG 
approximation and the red solid curve shows the result obtained with the LG+ approximation. 



Figure 6. Gluon-gluon scattering via gluon exchange 


This close agreement between the LC approximation and the LC-|- approximation 
seems to suggest that color flow issues are not important for the gap fraction. However, 
we believe that such a view would be mistaken and that both the LC approximation and 
the LC-|- approximation are inadequate for this problem. We can understand some of the 
issues without going into detail as follows. 

The rate of decrease of / as plp* decreases is controlled by the color flow in the event. 
To see this from the point of view of a dipole based parton shower, consider gluon-gluon 
scattering, as depicted in figure 6. Partons “a” and “b” scatter to produce final-state 
partons 1 and 2. Parton “a” has infinite positive rapidity, while parton “b” has infinite 
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Figure 7. A virtual correction to gluon-gluon scattering 

negative rapidity. In a leading color picture, partons “a” and “b” are color connected, as 
are “a” and 1, “b” and 2, and 1 and 2, as indicated in figure 6. There is another color 
configuration in which “a” is connected to 2 and “b” is connected to 1. 

These parton pairs form color dipoles. In a parton shower, the a-1 pair produces soft 
radiation in the angular region between pa and pi. That is, this radiation has large positive 
rapidity. Similarly, the b-2 dipole produces soft radiation with large negative rapidity. The 
a-b dipole produces soft radiation with any rapidity between a large positive value and a 
large negative value. Similarly, the 1-2 dipole produces soft radiation with any rapidity 
between the large positive rapidity of parton 1 and a large negative rapidity of parton 2. 
Thus it is very likely that the a-b, 1-2, a-2, and b-1 dipoles will produce radiation that fills 
in the gap. 

Now consider what happens in a parton shower based on the LC-|- approximation 
with perturbative corrections as in eq. (2.15). There can be a AV(r) contribution that 
represents a virtual gluon exchange between the two incoming partons, including imaginary 
contributions proportional to ivr. This gives a diagram in which a hard gluon is exchanged, 
then a soft gluon, as depicted in fig. 7. The exchange of two gluons is equivalent in color 
space to the exchange of an object with color 27, color 10, color 10, color 8 (in two ways), 
or color 1. With color singlet exchange, the two forward moving partons, “a” and 1, are 
completely disconnected in color space from the two backward moving partons, “b” and 2. 
If we get this color configuration in both the color ket state and the color bra state, we get 
forward radiation and backward radiation, but only a small amount of radiation into the 
gap region. Contributions like this can change the gap fraction. 

This same argument applies to virtual gluon exchange between partons 1 and 2. It 
also applies to virtual gluon exchange between “a” and 2 or between “b” and 1. It also 
applies to q-g scattering, q-q scattering, and q-q scattering. 

It is clear that this effect is not included in the LC-|- shower with no factors of AV(t). 
The existing literature [24-35] suggests, however, that the effect is important. To include 
it, we need to use the perturbative expansion of eq. (2.15). This expansion is not yet 
implemented in Deductor. 
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7 Conclusions and outlook 


We have examined the effect of color flow on the results from a parton shower genera¬ 
tor based on color dipoles, Deductor. On one hand, we used the leading color, LC, 
approximation that is standard for parton shower generators. We compared LC results 
to results from the LC-|- approximation that is available in Deductor. We used the 
LC-|- approximation up to terms suppressed by more than four powers of l/W- The LC-|- 
approximation has the advantage that it is exact in color for collinear emissions. It is, 
however, only approximate for wide angle soft gluon emissions. 

By comparing LC results to LC-|- results, we tested numerically how good the LC 
approximation is. We found that the LC approximation is quite good for the one jet 
inclusive cross section and for the distribution of the number of partons in a jet. There 
is an important caution needed here: one would really need to test the LC approximation 
against a full color treatment to be sure that the LC approximation is accurate for any 
given observable. However, testing the LC approximation against the LC-|- approximation 
certainly gives the LC approximation a chance to fail and it does not. We have tried some 
other observables and have found similar results. 

In section 3, we asked how much color suppression is typical. We found that color 
evolution beyond the leading color approximation is not rare. It is the norm. However, 
we find that simple observables are not much affected by deviations from leading color 
configurations. 

In one case, we found that the LC approximation agrees well with the LC+ approxi¬ 
mation, but we believe that neither approximation is adequate: the gap-between-jets cross 
section investigated in section 6. Our belief here is based on the analyses in the substantial 
current literature [24-35]: the color changing operators that need to be present are in fact 
not present in the LC-|- approximation. 

We believe that a parton shower approach to the gap-between-jets problem holds some 
promise because a parton shower approach can sometimes provide numerical results to com¬ 
plicated problems that are too unwieldy to solve analytically. Certainly, a parton shower 
can incorporate the effects of momentum conservation in the emission of soft gluons. These 
effects are not easy to incorporate in an analytical approach and they can be important 
[35]. In addition, a parton shower automatically includes potential radiation into the gap 
from a gluon that was radiated out of the gap. What is needed is a better treatment of 
color within the parton shower. In ref. [12], we described how one can insert perturba- 
tively the evolution operators that represent the difference between exact color and the 
LC-|- approximation. See eq. (2.15). It is perhaps not necessary to have a large number of 
insertions in order to have a reasonably good approximation: one virtual gluon exchange 
in the ket state and one in the bra state may be enough, as suggested by the numerical 
investigations in ref. [36]. We have not, however, implemented such a program, so we do 
not know how well it will work. 

Our tests do not include hadronization. We anticipate in the future linking Deductor 
to the string hadronization model [37, 38] of Pythia. For this purpose, we need to specify 
the probability for a given color density matrix element |{c}m)({c'}m| from eq. (2.1) to 
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correspond to a given classical string state. A method for doing this is given in section 8 of 
ref. [12]. We expect that the string configurations produced with the LC+ approximation 
may have more kinks than those produced with the LC approximation. Presumably this will 
necessitate retuning the parameters of the hadronization model compared to the Pythia 
default tune. It is difficult to predict what effects may remain after retuning. 

Finally, we note that a shower based on the LC+ approximation can start with any 
basis element |{c}m)({c^}m| of the color density matrix, including one with {c'}m /{c} m* 
This is not the case with the LC approximation as usually formulated. This is advantageous 
for matching to a perturbative NLO calculation since one can use the full NLO color- 
ordered amplitudes. There is a matching program [39] that matches to Deductor in 
an MC@NLO style. As presently implemented, this calculation reduces the NLO density 
matrix to {c'}m = {c}m, but we expect that this can be generalized to {c!}rn 7^ {c}m- 

Acknowledgments 

This work was supported in part by the United States Department of Energy and by 
the Helmholtz Alliance “Physics at the Terascale.” We thank Jeff Forshaw and Mrinal 
Dasgupta for helpful conversations about rapidity gap physics. 


A Using the leading color approximation 


Deductor uses the LC-I- approximation for color, as explained in some detail in ref. [12]. 
However, we may want to use just an LC shower. In particular, it is useful to turn the 
LC-|- shower off if the amount that the color suppression index has increased during the 
shower (after the hard interaction) reaches a preset value, (A/)max- After that, we can 
proceed with an LC shower. The leading color approximation is commonly used in parton 
shower algorithms, so it would hardly need an explanation. However, we need to run an 
LC shower starting with a color density matrix 


^ p{{p,f,c',c}m,t)\{c}m){{c'}m\, (A.l) 


where possibly {c'}m 7 ^ {c}m- This is straightforward, but it is not commonly done, so 
we explain it in this appendix. In order to keep the appendix brief, we assume that the 
material in refs. [8] and [12] is known. 

The method for switching from LC-|- to LC that we explain here and that is imple¬ 
mented in Deductor is different from the method outlined in appendix C of ref. [12]. 

The LC approximation is based on replacing the color group SU{Nc) by U{Nc). There 
are gluons, labeled with q = {1, 2,..., N^}. The generator matrices t^,..., in 

the fundamental representation are the usual ones. The new generator is 


f ^ — 




(A.2) 


- 15 - 




This gives the identity 


= l ■ (A.3) 

a=l 

That is, aside from the factor 1/2, we can think of each gluon as carrying one color Nc line 
and one color Nc line. The factor Tr = 1/2 comes from our conventional normalization of 
the generator matrices, 

Tr[t“t'] = lSab. (A.4) 

With the dehnition Eq. (A.2), this normalization remains true for all a and b. 

We define the structure constants fabc that give the couplings of the gluons to each 
other by 

= i/afect'" ■ (A.5) 

With this definition, we see that the U{1) gluon that we have added, with index Nc, does 
not couple to the other gluons. That is, fabc = 0 if any of a, b, or c is N^. 

Our analysis makes use of the standard color basis dehned by two sorts of vectors, as 
in ref. [8]. First, there are open string vectors 

= iV(S)-^/2[^a2^a3 ... . (A.6) 

Here oi is a quark color index, an is an antiquark color index, and the other ai are gluon 
color indices. The N are U{Nc) color matrices in the fundamental representation. Also, 
N{S) = Nc{Nc/2)^~‘^ is a normalization factor such that = 1. Second, there are 

closed string vectors 

^(0)W = Tr[t“H“2 • • • . (A.7) 

Here all of the a, are gluon color indices. Again, N(S) = {Nc/2)'^ is a normalization factor 
such that ('I'I'k) = 1. (With SU{Nc), is slightly different from 1.) The most general 

color basis vector, which we denote by |{c}m), is a product of these two kinds of units. 

The color basis states are normalized to ({c}m|{c}m) = 1- They are not, how¬ 
ever, generally orthogonal. However, when {c}m and {c'}m are different, one hnds that 
({c'}m|{c}m) = C>(l/-Ac) or 0(1/Nc). That is, the basis vectors are orthogonal only in 
the Ac —> oo limit. 

The evolution equation for the shower state as a function of the shower time t has the 
form of a linear equation 

4|p(i)) = !«,(*)-V(()]|p(*)) . (A.8) 

Here 'Hi{t) is the parton splitting operator. It is given for full color in eq. (5.7) of ref. [12]. 
The LC-t- approximation is defined in ref. [12] by restricting the possibilities for color states 
that can be reached by a splitting. We can define the LC approximation by restricting these 
choices further. This is illustrated in figures 8 and 9. 

In hgure 8 we illustrate the color structure for the splitting of a gluon into two gluons, 
numbered 1 and 3. The same splitting occurs in both the bra and ket states. Note that in 
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Figure 8 . Identity for color dependence of splitting of gluon 1 in both the bra state and the ket 
state in the case {c'}m 7 ^ {c}m- In the LC+ approximation, one keeps all four terms. In the LC 
approximation, we keep only the first two terms. 



Figure 9. LC+ approximation for the splitting of gluon 1 in the ket state with the participation 
of helper parton 2 in the bra state in a case with {c'}m 7 ^ {c}m- In the LC approximation, we keep 
only the first term. 


this example, the color structures of the bra and ket states are not the same. Triple gluon 
vertices are ifabc- Using eq. (A. 5), we can expand the left hand side of figure 8 into a linear 
combination of four products of color basis functions, as shown in the right hand side of 
figure 8. In the LC+ approximation, we keep all four of these terms. We define the LC 
approximation to keep just the first two terms. In the shower Monte Carlo algorithm, we 
choose one state or the other with probability 1/2. 

In figure 9 we illustrate the color structure for the radiation of a gluon 3 from gluon 

1 in the ket state, but with interference with the radiation of gluon 3 from “helper gluon” 

2 in the bra state. Expanding the left hand side into products of basis functions produces 
four terms. Of these, the LC+ approximation keeps the two terms shown in the right hand 
side of figure 9. We define the LC approximation to keep only the first term. 

The general principle for emission of a gluon should be clear from these two examples: 
the new gluon (3) must not cross the emitting gluon (1) in the color diagram. For emission 
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of a gluon from a quark, the situation is simpler: there is only one allowed structure in the 
LC+ approximation (as in figure 16 of ref. [12]) and this structure is retained in the LC 
approximation. There is also the possibility of g —>• g + g splittings. With the color group 
U{Nc), there is only one possible new color state after the splitting. The new state is given 
by connecting the new quark and antiquark to the previous color string by using eq. (A. 3). 

We also need the operator V{t) in eq. (A. 8). This operator leaves the number of 
partons unchanged. Its exponential is the Sudakov factor that appears between splittings 
in the shower. We dehne V{t) from %i{t) in such a way that the inclusive cross section 
that we started with at the hard interaction is not changed at all by the shower. The 
operator V{t) contains a color factor denoted by N{k,l,{f}m+i) in ref. [12]. For the 
LC+ approximation, it is given in eq. (6.12) of ref. [12]. For the LC approximation, we 
need to calculate N{k,l,{f}m+i) using the color group U{Nc) and including only the LC 
contributions. For all cases of gluon emission except for a g —g + g splitting as in figure 
8, N{k,l, {f}m+i) is Nc/2, where the W comes from having one more “quark” loop and 
the 1/2 is the 1/2 in eq. (A. 3). For a g —)• g + g splitting, there are two terms, so the net 
color factor in V{t) is 2{Nc/2) = W = ^A- For a g — >■ g + g splitting, N{k,l,{f}m+i) is 
just Tr = 1/2. 

For the case of a g —)• g + g splitting in either the ket state or the bra state or both, the 
color factor in V{t) is Nc/2 when we calculate this way. One normally uses Cp instead. We 
arrange to obtain the conventional net factor of Cp instead of Nc/2 by simply multiplying 
the splitting probability by 2Cp/Nc in for the LC approximation. 

We need one more adjustment. When the change in color suppression index during the 
shower reaches (A/jmax) we switch from the LC+ approximation to the LC approximation 
as described above. This means changing the rules for calculating the color factor in quan¬ 
tum probabilities. If the color state before this switch is |{c}m)({c'}m| then the correspond¬ 
ing contribution to the inclusive cross section contains the color factor {W}m\{c}m)y 
where the subscript indicates a calculation using the color group SU{Nc). After the switch, 
we have the same nominal color state, |{c}m)({c^}m|, but we have changed the definitions 
so that now the inclusive cross section contains the color factor ({c'}m j calcu¬ 

lated using the color group U{Nc). To preserve probabilities, we multiply the weight for 
the current shower state by 


Cn = 




(A.9) 


This factor is normally close to 1, but in some situations the numerator can be zero. We 
note that Deductor contains fast algorithms for calculating both the numerator and the 
denominator. 


B Another approach 

Platzer and Sjodahl have advocated another approach [13] to the problem of incorporating 
quantum color in a parton shower. This approach is dramatically different from that of 
refs. [8, 12]. 
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Recall that for a parton shower with quantum color, the partonic state \p{t)) is a 
density matrix in the quantum color states. How ought \p{t)) to evolve in a parton shower 
with full color? In our opinion, the most basic requirement is the following. At the end 
of the shower at shower time tf, we have a partonic state \p{t{)) that can be expanded 
in powers of Os- The first order term in this expansion should be given by the Born 
approximation, |/5 b)) to |/9) convoluted with collinear/soft approximations to one loop 
QCD graphs with an infrared cutoff corresponding to the final shower time tf. Ref. [13] 
works with \p{t)^, but the first order term in the expansion of \p{ti)) is not given by |pb) 
convoluted with collinear/soft approximations to one loop QCD graphs with their correct 
color structure. 

We can be more specific. Following the structure of quantum statistical mechanics and 
the idea that in a parton shower one parton can split into two partons, the evolution of the 
shower state can be represented as a linear equation, as in eqs. (2.7) and (2.8) [8]. This 
equation involves an operator that increases the number of partons by one and an 

operator V{t) that leaves the number of partons and their momenta and flavors unchanged. 

Ref. [13] has T-Li{t) in agreement with ref. [8]. The color structure of V{t) is very well 
known [8, 24-34, 36, 40]. However, in ref. [13] with “full color,” the operator V{t) is absent.^ 
It is replaced by a Sudakov factor that does not change the partonic color state. The “full 
color” treatment of ref. [13] is perhaps an approximation to a full color treatment, but it 
is not clear to us what sort of approximation it is or whether the approximation can be 
improved. 
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